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Abstract. In this paper we present a quantitative analysis of the chemical abundances in "monolithic" N-body- 
Tree-SPH models of elliptical galaxies to determine whether the ratio [a/Fe] increases with the galaxy mass as 
suggested by observational data. 
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1. Introduction 

Twenty years after the seminal papers on the formation 



of elliptical galaxies (EGs) by Larson (1975) and Toomre 



(1977) there is still much disagreement among astronomi- 
cal community on both the process of formation and evo- 
lution of early- type galaxies. 

Are EGs formed in a single collapse during a short 
epoch of activity (monolithic scheme) or continuously by 
hierarchical mergers at similar levels by several mergers 
(hierarchical scheme)? Once created, is the stellar popu- 
lation of EGs just passively evolving or do minor merg- 
ing/accretion events drastically change their characteris- 
tics frequently? What is the influence of density environ- 
ment? 

In the monolithic scenario EGs are supposed to suffer 
dominant star forming activity and consequent chemical 
enrichment at very early epochs followed by quiescent evo- 
lution, marginal episodes of star formation of both inter- 
nal and external origin may occur under suitable circum- 
stances (e.g. interactions). 

In the hierarchical scenario EGs are supposed to be 
formed by mergers (one to several) of smaller units, each 
episode inducing star formation and chemical enrichment. 

Both scenarios are able to reproduce only part of 
the observed properties of EGs. For instance the "mono- 
lithic scheme" , on which current semi-analytical models of 
chemical and spectro-photometric evolution of EGs (see 
below) are based, cannot account for the wide morpho- 
logical and kinematical variety of EGs, for instance the 
presence of counter-rotating cores, small disks, recent stel- 



lar activity ( 


Schweizer 1990|; pchweizer & Seitzer 1992 


Longhetti et al. 1998a, 




2ooq). 



In contrast, despite its success in modeling the dy- 
namical structure of EGs, the "hierarchical scenario" suf- 
fers a point of difficulty as far as the chemical properties 
are concerned. Indeed, EGs often have different ratios be- 
tween the abundances of some chemical species. They are 
in general more metal-rich and enhanced in a-elements 
with respect to spiral galaxies. In the merging scenario 
this could be justified by an additional chemical enrich- 
ment with strong enhancement in a-clcmcnts must take 
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place. See the reviews by Chiosi (2000| , |200lD . 

Recently Chiosi & Carraro (2001) have investigated by 
means of N-Body, Tree-SPH models whether the mono- 
lithic scheme of galaxy formation and evolution may in- 
terpret the observational properties of EGs. They have 
studied in particular how the history of star formation de- 
pends on the initial conditions (density) and/or mass of 
the proto-galaxy showing that the monolithic scheme can 
indeed reproduce many fundamental properties of EGs, 
such as colors, mass to light ratios, the color magnitude 
relation (CMR), and the so-called Fundamental Plane (see 
[Chiosi fc Carraro 2001 , for all other details). 

In this study we intend to examine whether the same 
models would also predict chemical properties compat- 
ible with the observational data. In particular we esti- 
mate their degree of enhancement in a-elements and check 
whether they can account for the [a/Fe] -mass relation of 
EGs (see below). 

The plan of the paper is as follows. In Section 1 we 
discuss the internal contradiction between the kind of con- 
straint imposed by the CMR and the [a/Fe]-mass relation- 
ship. In Section 3 we present the main results from the dy- 
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namical models, outline the companion chemical models, 
and report on the results for the chemical abundances and 
abundance ratios. In Section 4 we present an ad hoc exper- 
iment aimed at pinning down the main parameter driving 
the level of a-enhanced in model galaxies. In Section 5 
we summarize a few basic relations defining the enhance- 
ment factor of a single stellar population (SSP) of assigned 
chemical composition [X,Y,Z], present a simple algorithm 
to calculate the spectral indices Hp, Mg2, (Fe) and Mgb, 
apply it to our model galaxies, and derive their mean in- 
dices. Finally, these latter are compared with the obser- 
vational data paying particular attention to the Mg2 vs a 
(velocity dispersion) relation and the vs Mg2 and 
vs [MgFe] planes. In Section 6 we draw some concluding 
remarks. 

2. A long lasting contradiction 

EGs show important properties that are used as key con- 
straints to theoretical models aimed at understanding the 
formation mechanism and the past evolutionary history. 
Chief among these are the CMR and the enhancement in 
a-elements suspected to exist in the brightest (most mas- 
sive) ellipticals. 

2.1. The Color-Magnitude Relation 

The colors of EGs get redder at increasing luminosity 



and/or mass ( 


Matthews & Baker 1971 




Larson 1974; 


Bower et al. 1992a 


b). The relation is tight for cluster 



for nearby field EGs and those in small groups ( Bchweizer 
fc Seitzer 1992|) . 



The CMR is con ventionally interpre t ed as a mass 

Dressier 



metallicity sequence ( |Fabcr et al. 1977 



1984 



Vader 1986): massive galaxies reach higher mean metallic- 
ities than less massive ones. Starting from the original sug- 
gestion by Larson (1974 ), the mass-metallicity sequence is 
attributed to the so-called supernova driven galactic wind 
(SDGW) mechanism halting star formation at an age that 
depends on the galaxy mass. In short, the overall duration 
of the star forming activity should be longer at increasing 
mass of the EGs (Atsp oc Mg)- 

The tightness of the CMR is thought of to reflect the 
age of a galaxy stellar population. The CMR of the Virgo 
and Coma clusters, where the dispersion is very small, is 
compatible with most of the stars in EGs being formed at 
redshift z > 2 which roughly corresponds to a lookback 
time of about 10 Gyr for a Universe with go = 0.5 and 
Hq = 70 km s^^ Mpc~^. The absolute value of the age is 
less of a problem here. Nevertheless, secondary bursts of 
star formation at z = 0.5, forming about 10% of the stellar 
mass could be also accommodated ( Bower et al. 1992a| ,|b|). 
T he possibility of secondary ac tivity has also been invo ked 
by Kuntschner fc Davies (199g ) and Kuntschncr (2006|) on 



a different ground. In contrast the more dispersed CMR of 
the field and small group EGs is compatible with ages of 
the component stars spread over several billion years. This 



has also been suggested using the line index diagnostic 
by Gonzalez (1992) and Trager et al. (2000a ^ with a 
tendency for lower velocity dispersion (mass) systems to 
be younger (see also 



Bressan et al. 1996| ). 



2.2. The [a/Fe]-Mass Relation 

Gradients in line strength indices Mg2 and (Fe) (and oth- 
ers) have been observed in EGs ( Worthey 1992|; ponzalez 
1993|; pavies et al 199"3| paroUo et al. 1993|; paroUo fc 



1995 



Danzigcr 1994a| ,|b|; |Balceils fc Peletier 1994| ; [Fisher et al 



199£). The hint arises that the indices are stronger 
toward the center. There is also the controversial ques- 
tion whether t he gradients in Mg2 and (Fe) are the same. 
Accor ding to ( CaroUo et al. 1993 ; CaroUo fc Danzigei 
1994a| ,|b|) the gradient in Mg2 is steeper than the gradi- 
ent in (Fe), which is interpreted as indicating that the 
abundances of a-elements (Mg, O, etc.) with respect to 
iron are enhanced - [a/Fe] >0 - in the central regions. 
However, Kuntschner (1998|) reaches the opposite conclu- 



sion and pavies et al. (200l| ) show for NGC4365 that the 



ratio [Mg/Fe] stays constant within the data range. 

Furthermore, limited to the nuclear regions, indices 



vary passing from one galaxy to another (Gonzalez 1993; 
Trager et al. 20004 |b|). Looking at the correlation be- 
tween Mg2 and (Fe) (or similar indices) for the galaxies in 
the above quoted samples, Mg2 increases faster that (Fe) 
which is once more interpreted as due to enhancement of 
a-elements in some galaxies. In addition to this, since the 
classical paper by Burstein et al. (1988| ) , the index Mg2 is 
known to increase with the velocity dispersion (and hence 
mass and luminosity) of the galaxy. Standing on this body 
of data the conviction arose that the degree of enhance- 
ment in a-elements ought to increase passing from dwarf 
to massive EGs (|Faber et al. 1992|; [Worthey et al. 1994 ; 
Matteucci 1994], |1997|; fVIatteucci et al. 1998|). 



Understanding these properties rests on stellar nucle- 
osynthesis and star formation. In brief, all heavy elements 
observed in the interstellar medium (ISM) are produced 
and ejected by stars, in particular via supernova explo- 
sions. Furthermore, iron is mainly produced by the long- 
lived Type la supernovae (accreting white dwarfs in binary 
systems in the most popular scheme) whereas a-elements 
are mainly generated by the short-lived Type II super- 
novae. Therefore, to enhance the a-elements there are sev- 
eral avenues: 

— (i) A different ratio between Type la and Type II su- 
pernovae. 

— (ii) A different time scale of star formation. 

— (iii) A prc-enriched interstellar medium. 

The simplest and most widely accepted interpretation 
is the one based on the different duration of the star form- 
ing period and the different contribution to chemical en- 
richment by Type la and Type II supernovae. In this con- 
text, it is worth recalling that the minimum mean time 
scale for a mass accreting white dwarf in a binary system 
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to get the supernova stage is >0.5 Gyr. Therefore, the iron 
contamination by Type la supernovae occurs later as com- 
pared to the ones of Type II dGreggio fc Renzini 198^ ). 
It follows from this that with the standard SDGW model 
and classical initial mass function (IMF), the time scale 
of star formation and galactic wind must be shorter than 
about 1 Gyr not to decrease the initial [a/Fe]>0 (when a- 
elements are mostly produced) to [a/Fe]<0 (when iron is 
predominantly ejected). In other words, to reproduce the 
observed trend of the [a/FeJ-mass relationship, the total 
duration of the star forming activity ought to scale with 
the galaxy mass according to AtsF oc M^^, which opposes 
to the expectation from the CMR. 

This point of contradiction cannot be avoided by the 
standard chemo-sp ectro-photometric mo dels - both the 
closed-box ones by Bressan et al. (1994 ) a nd those with 



infall of primordial gas at early epochs by Tantalo et al. 
(1996171 |l 998"!) - because they aU stand on the SDWG 



scheme and are tailored to reproduce the CMR. Therefore, 
they fail to fit the companion mass-a-elements relation- 



ship. This issue has been addressed by Tantalo fc Chiosi 
(200l|)]who have carefully examined the effect of key pa- 
rameters of the classical models, namely the type and ef- 
ficiency of star formation and the role played by the rate 
of infall on the degree of enhancement in a-elements. The 
main result of this analysis is that the duration of star 
formation activity and age at which galactic winds set in 
do not bear very much on the degree of a-enhancement. 
This in contrast is determined by the detailed time depen- 
dence and intrinsic efficiency of star formation. In other 
words, a galaxy whose stellar activity predominately oc- 
curs at very early epochs (say within the first Gyr), even 
if followed by a long tail of star formation at smaller ef- 
ficiency (it does not matter whether continuously or in 
recurrent episodes), would nowadays appear as enhanced 
in a-elements. Obviously, the opposite would hold for a 
galaxy whose star formation is diluted at low, yet compa- 
rable levels of intensity over long periods of time. 

Finally, another short comment is worth here. 



Following Ferrini (2000 ), the above point of contradiction 
could stem from the hypotheses that (i) SFR and con- 
sequent rate of supernova explosions (mostly type II) are 
proportional to the gas content (and total galaxy mass AIg 
in turn); (ii) the binding energy of the gas is proportional 
to Mg; and finally (iii) the star forming activity stops after 
the galactic wind has occurred. All of the above hypothe- 
ses are oversimplified descriptions of reality, for instance 
SF is known to occur in molecular clouds, which are a sub- 
component of the gas content (in other words the SFR can 
be a non linear function of the gas content), SF not neces- 
sarily halts after the first galactic wind episode and finally 
the galactic winds can spread over significant intervals of 
time. All this could destroy the simple relation between 
the duration of SF and galaxy mass and perhaps allevi- 
ate the di screpancy be tween the CMR and the a-element 
problem ( Ferrini 2000| ). 

The present approach based upon realistic dynamical 
models of early type galaxies, in which neither star forma- 



tion comes to an abrupt halt nor galactic winds suddenly 
set up, see Chiosi & Carraro (2001) for all details, is per- 
haps a first step toward a self-consistent formulation of 
the problem. 

3. Enhancement of a-elements in dynamical 
models 

3.1. The dynamical models 



As we have already anticipated, Chiosi fc Carraro (2001 ) 
adopting the monolithic scheme have calculated fully hy- 
drodynamical models of early-type galaxies with different 
total mass Mt. The models included dark and baryonic 
mass, Md and Mb respectively, in cosmological propor- 
tions (Md = 0.9 X Mt and Mb = 0.1 x Mt), star forma- 
tion, chemical evolution, energy feed-back, and cooling by 
radiative processes. 

In order to analyze the effect of the initial density of 
the proto-galaxy three sets of models were considered: (i) 
a first group (labeled A) whose initial density is given 
by (p)o — 200 X Pu{z) (where /9«(z) is the mean density 
of the Universe as a function of red-shift); (ii) a second 
group (labeled B), in which the initial density has been 
arbitrarily set equal to {p)o ~ Pu{z); (hi) finally a third 
group limited to a unique value of the total mass (Mt = 
10^ M0) in which the initial density has slightly varied 
above and below the case B value. They are hereinafter 
indicated by LD (lower density) and HD (higher density). 
For all other details the reader to refer to phiosi fc Carraro 
|(200lD . The main results of this study can be summarized 
as follows: 

(i) Independently of the total mass, galaxies of high ini- 
tial density (case A), undergo a prominent initial episode 
of star formation ever since followed by quiescence. 

(ii) The same applies to high mass galaxies of low ini- 
tial density (case B), whereas the low mass ones undergo 
a series of burst-like episodes that may stretch over a con- 
siderable fraction of their lifetime. The details of their star 
formation history are very sensitive to the value of the ini- 
tial density as shown by the models of the third group. 

(iii) The mean and maximum metallicity increase with 
the galaxy mass. Therefore these models can account for 
the CMR of EGs. 

(iv) The model galaxies suffer galactic winds ejecting 
conspicuous amounts of mildly enriched gas into the intra- 
cluster medium. However, they do not follow the simple 
SDGW mechanism envisag ed by ^ arson (1974| ): when the 
total thermal energy of the gas heated up by supernova 
explosions (radiative cooling included) exceed the gravi- 
tational potential energy of it, the remaining gas is sud- 
denly ejected and star formation is halted (see for instance 
Bressan et al. 1994; Tantalo et al. 199G, 1998 for more de- 



tails on this type of models). In dynamical models the sit- 
uation is more complicated. Looking at radial velocity Vr 
of the gas and the escape velocity Vosc, at any time there 
is a layer a which some gas particles may meet the condi- 
tion Vr > Vcsc- All these gas particle can freely leave the 
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Table 1. Input parameters of the Tree-SPH models of EGs by Chiosi & Carraro (2001). 



Model 


Mt 


(P)0 


ATsF 


TpSF 


(*> 


M, 


0"B 








MQ/kpc^ 


Gyr 


Gyr 


Mo/yr 


Mq 


km s^^ 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 






Models with {p)o — 


200 X pu{z) 








lA 


1 X 10^^ 


4.2 X 10^ 


0.5 


0.25 


3200 


6.4 X 10" 


530 


99 


6A 


1 X 10'' 


4.2 X 10^ 


0.3 


0.15 


0.07 


3.2 X 10^ 


15 


99 








Models with (p)o 


— Pu{z) 








IB 


5 X 10^^ 


3.8 X 10'' 


2.5 


2.0 


1760 


4.4 X 10^^ 


464 


99 


2B 


5 X 10^^ 


3.8 X 10" 


2.5 


1.8 


160 


4.0 X lO" 


210 


98 


3B 


1 X 10^^ 


3.8 X 10-* 


2.5 


3.0 


25 


6.4 X 10^" 


97 


98 


5B 


1 X 10^0 


3.8 X 10-* 


3.5 


2.7 


0.11 


4.0 X 10** 


19 


94 


6B-LD 


1 X 10'' 


5.5 X 10^ 


6.5 


4.0 


0.002 


1.3 X 10^ 


3 


44 


6B 


1 X lO'' 


3.8 X 10'' 


4.5 


2.5 


0.005 


2.5 X 10^ 


7 


76 


6B-HD 


1 X lO'' 


5.8 X 10" 


3.5 


1.9 


0.022 


7.7 X 10^ 


5 


91 


7B 


1 X 10* 


3.8 X 10" 


6.5 


3.0 


0.0003 


1.8 X 10'' 


4 


90 



galaxy in form of galactic wind. Therefore galactic winds 
are a kind of process taking place continuously, at least 
as long as there are gas particles reaching and exceeding 
the above condition. Furthermore, the total amount of gas 
lost in the wind AMgw depends on Mt- The following use- 
ful relations are found AMgw/Me = -0.091ogMT + 1.256 
for models B and AMgw/Me = -O.OeiogMx + 0.99 for 
models A. Low mass galaxies are somewhat more efficient 
emitters of partially processed material. 

(v) How does it happen that low mass galaxies may 
eject large quantities of gas and yet form stars for long 
periods of time? In a normal or giant elliptical whose 
gravitational potential well is deep, the gas tends to re- 
main inside the potential well despite being heated up by 
energy injection by the supernova explosions and stellar 
winds accompanying star formation. A sort of equilibrium 
is reached in which gas heating and cooling by radiative 
processes balance each other, forcing gas to consumption 
by star formation. In contrast, in a low mass galaxy, ow- 
ing to the shallower gravitational potential well, the same 
processes heat up gas, part of it spills over the poten- 
tial well and part flies to very low density regions so that 
star formation in the inner regions declines. Eventually, 
the gas pushed away to large distances but still bound to 
the galaxy cools down and falls back into the central fur- 
nace thus re-fueling star formation. Therefore, this latter 
keeps going for long periods of time in a series of burst- 
like episodes till either gas is exhausted or what remains 
is heated up so much (its gravitational binding energy in 
meantime gets smaller and smaller compared to the ther- 
mal energy) that it eventually leave the galaxy. 

(vi) Finally, thanks to their star formation history 
these models have the potential of improving upon one 
of the most embarrassing difficulties encountered by the 



standard supernova driven galactic wind (SDGW) model 
as far as the chemical abundances and abundance ratios 
are concerned, i.e. the high [a/Fe] ratio inferred in mas- 
sive EGs and the long duration of the star forming period 
required by the SDGW model to explain the CMR. 

In Table [l| we summarize a f ew properties of the mod- 
els by [Chiosi fc Carraro (2001 ) that are relevant to this 
study (we follow the sa me identification coding used by 
[Chiosi fc Carraro (2001 ) for the sake of an easy compar- 
ison). Column (1) identifies the model. Column (2) gives 
the total mass Mt in solar units. Column (3) gives the 
mean initial density (p)o in units of Mq kpc~^. Columns 
(4) and (5) list the total duration of the star formation ac- 
tivity ATgF in Gyr, and the age Tpgp in Gyr at which the 
maximum stellar activity occurs, respectively. Column (6) 
lists the mean star formation rate (^) in units of Mq/yt. 
Column (7) gives the present-day mass in stars in solar 
units. Column (8) contains the central velocity dispersion, 
(7 of baryonic matter in kms~^. 



3.2. The detailed chemistry 

Although the description of the chemical enrichment 
adopted by [Chiosi fc Carraro (2001 ) was fully adequate 
to the purposes of N-body Tree-SPH simulations, the de- 
tails of chemical abundances and abundance ratios were 
out of reach owing to the numerical complexity of those 
calculations. 



To copy with this obvious limitation of the dynamical 
simulations, we adopt here an indirect approach which 
albeit simple is sufficiently accurate, i.e. we insert the SFR 
of dynamical models into the chemical-evolution code of 
Portinari fc Chiosi (1999|) . 
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Since chemical evolution affects only the baryonic mass 
of the galaxy it is sufficient to follow the evolution of an 
ideal galaxy having the same baryonic mass Mb, initially 
in form of gas. The chemical models are in the so-called 
cZoserf-6oa; formulation (see Tinsley 198C| ) as the dynamical 
ones. Finally, care has been paid to secure that the yield 
of heavy elements and gas per stellar generation are the 
same (or nearly the same) as those used in the dynamical 
models so that the present-day metallicity and mass in 
stars of the dynamical and chemical models with the same 
baryonic mass and star formation history coincide. 

Sketch of the chemical models. The equation governing 
the abundances of elemental species in chemical models 
of galactic evolution are well known so that no details are 
given here. The reader is referred to Portinari et al. (1998 ) 
for the detailed description of the equations used also in 
the present study. Suffice here to recall a few important 
ingredients: 

(i) The oquationo arc caot ao function of the dimonoion 



The rate of star formation. The rate of star forma- 
tion (SFR) of the dynamical models is according to the 
^chmidt (1959 ) law applied to each gas particle. Therefore, 
the total star formation rate is also the Schmidt law ap- 
plied to the whole system: 



dMg(t) 
dt 



^Mg(t)« 



(1) 



where v—1 and k—1 as in Chiosi fc Carraro (2001] ). The 
SFR of the dynamical models must be rescaled accord- 
ing to the normalization adopted in the equations of the 
chemical models. Upon normalization, the star formation 
rate becomes: 



-r , ^ d , M„ , 



(2) 



less va riabloD 



The initial ma ss function. We have ad opted he re the 
same IMF used by [Chiosi fc Carraro (200ll ), i.e. the |Millei 
fc Scalo (197S| ) over the mass interval Mi=0.1 to Mu=120 



G,(t).MiW, G.(t) = M^. 

' Mb ^ ^ Mb 



Gg,i(t) = Gg(t) X Xi(t) 



where Xi(t) are the abundances by mass of the elemental 
species i^ By definition EiXi(t)=l. 

(ii) The restitution fractions of the elements i from a 
star of mass M at the end of its life are obtained by means 
the matrix formalism of (see Talbot & Arnett 1971, 1973| , 
1975| ). 



Mq . It is worth of interest to look at the fractionary mass 
stored by the IMF in the mass interval whose stars ef- 
fectively contribute to nucleosynthesis in a time interval 
shorter than the age of the galaxy. A suitable choice is the 
mass range Mn^lM© to Mu. The above fraction ( is 



/^;^(M)xMxdM 



/^"0(M)xMxdM 



(3) 



(iii) No use of the instantaneous recycling is made (i.e. 
we take the life time tm of stars of mass M into account. 

(iv) We consider the separated contributions to nucle- 
osynthetic ejecta by single and binary stars. Single stars 
can terminate their life quietly becoming white dwarfs and 
ejecting part of their mass in form of stellar winds or ex- 
ploding as Type II supernovae. Accreting white dwarfs in 
binary stars are thought of to b e the site of Type la super- 
nov ae as long ago suggested by Greggio fc Renzini (1983 ) 



and Matteucci fc Greggio (1986 ). To evaluate the number 
of Type la supernovae we need do fix the minimum (MBm) 
and the maximum (MBj\/) value for the total mass (Mbm) 
of the binary systems, the distribution function (f(/i)) of 
their mass ratios, where fj-mm is the minimum value of 
it, and finally the fraction A of binaries with respect to 
the total. It is assumed here that binary stars as a whole 
obey the same IMF of single stars. We adopt B,„ — SM©, 
Bm = 12M0, and A = 0.02. 



(v) The stellar ejecta are from Marigo et al. (1996 
1998| ) and Portinari et al. (1998), and incorporate the de- 



pendence on the initial chemical composition as described 



by Portinari et al. (1998) 



(vi) Finally, the stellar lifetimes tm are from Bertelli 



et al. ( L994). We take also into account the dependence of 
Tm on the initial chemical composition. 

^ Fifteen chemical elements are followed in detail, namely H, 
*He, i^c, "c, i^N, i^O, 20Ne, ^^Mg, ^^Si, =^^8, «Ca, ^^Fe, and 
the isotopic neutron-rich elements nr obtained by a-capture on 
"N, i.e. i«0, ^^Ne, ^^Mg. 



With the slope of the [Miller fc Scalo (1979D IMF and the 
above mass limits, ( 0.25. 

3.3. Chemical abundances and abundance ratios 

The chemical evolution of all the models presented in 
Table |^ has been followed and the results of interest here 
are best represented by the series of Figs. |l| through ^. 

In Fig. we show the time variation of the SFR 
(top panel), the metallicity Z (mid panel), and the abun- 
dance ratio [0/Fe] (bottom panel) for the models 3B with 
Mt = lO^M© and 6B-LD with Mt = IO^Mq as indi- 
cated. As expected, during the period of stellar activity, 
the metallicity increases, whereas when star formation is 
over the metallicity remains either constant or slightly de- 
creases by dilution with gas ejected by old stars of low 
metal content. Looking at the time variation of the abun- 
dance ratio [O /Fe] ~ a measure of the a-enhancement level 
- we note that 

(i) Within the first 0.1 Gyr the abundance ratios are low- 
ered to values comparable to the mean ones under the 
contamination by the first Type la supernovae. 

(ii) Even a small amount of stellar activity changes the 
initial abundance ratios to much lower values. 

(iii) Subsequent star formation does not alter the abun- 
dance ratios significantly. See for instance the case of the 
10^^ Mq galaxy in which the prominent burst occurring 
from 2 to 3.5 Gyr changes (increases) the [0/Fe] ratio by 
less than A[O/Fe]=0.2 (less than a factor of two). 
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Fig. 1. Bottom: the temporal variation of the abundance 
ratio [0/Fe] in models 3B and 6B-LD. Middle: the tem- 
poral variation of the metallicity in the same models. Top: 
the SFR of the same models normalized to the peak value 

SFRMax- 



(iv) Once star formation has ceased, almost immediately 
the effect of Type II supernovae stops and that by Type la 
takes over first rapidly and later slowly decreasing [0/Fe] 
towards the solar value or below this. 

In order to answer the question whether or not the 
present-day stellar content of these galaxies would appear 
to be enhanced in a-elements, we look at the cumulative 
distribution of the number of living stars in different bins 
of [0/Fe]. From an operational point of view a galaxy will 
be considered dominated by an a-enhanced stellar content 
if X;N([0/Fe] > 0)/Nt > 0.5, where EN([0/Fe]) is the 
number of stars in bins with [O/Fe]>0 and Nt is their 
total number. The opposite if N( [0/Fe] < 0)/Nt < 0.5. 
The situation for the two galaxies under examination is 
shown in Fig. | (lO^^ Mq) and Fig. | (10^ Mq). 

The massive galaxy (Fig. ||), characterized by a rather 
short (~ 2.5 Gyr) and prominent burst of star formation 
activity has about 98% of its stellar content in bins with 
with [0/Fe]> 0. Whereas the low mass, low initial den- 
sity galaxy (Fig. ^), whose star formation stretched over 
long periods of time in several recurrent episodes, has the 
majority of stars in bins with [0/Fe]< (~44%). 

The percentage of stars with N([0/Fe])/NT > 0.5 are 
listed in column (9) of Table |l|. Finally, in columns (3) 
and (4) and of Table ^ to better described below we list 
the maximum, ZMax, and the mean metallicity, (Z), of our 
model galaxies. 



Fig. 2. The top panel shows the SF history of the 3B 
model with total mass Mt = IO^^Mq and baryonic mass 
Mb = 10"Mq. The SFR is in Mg/yr and time is in Gyr. 
The bottom panel gives the cumulative number of living 
stars per bin of metallicity Z (so/id line and bottom x-axis) 
and per bin of [0/Fe] {dotted line and top X-axis). The 
shaded zone shows the region populated by the stars with 
[O/Fe]>0, the a-enhanced component. The percentage of 
star numbers with [O/Fe]>0 is about 98%, whereas that 
with [O/Fe]<0 is 2 %. 



4. What does determine the enhancement of 
a-elements? 

The entries of Table |l| show that the percentage of stars 
with [O/Fe]>0 decreases with the total duration of the 
star formation activity ATsf- However all models but 
for the case 6B-LD have percentages of a-enhanced stars 
higher than 50%. In order to understand under which cir- 
cumstances (shape of the SFR and ATsf) a galaxy would 
build a stellar population with a percentage of a-enhanced 
stars less than 50%, we have performed the following ex- 
periment. We have taken one of the most unfavorable 
cases, i.e. a typical galaxy with Mt = 10^^ M© and hence 
Mb = 10^^ Mq, imposed that the same amount of mass 
in stars is build up at the end of the star forming activ- 
ity, and varied both the shape and duration, ATgp of the 
SFR. 

We adopt a simple analytical expression for the SFR 
as illustrated in Fig. ^. This is idealized to a trapezium 
of sides A, B and ATsf- By varying the ratio B/A, the 
shape of the SFR changes from case (1) in Fig. ^ with 
B/A<< 1 (SFR peaked in the past), to a constant SFR 
with B/A=l of case (2) in Fig. |, and finally to case (3) in 
Fig. ^ (SFR peaked toward the present). The total area of 
the trapezium (rectangle) yields the total mass Mg in stars 



R. Tantalo & C. Chiosi: Enhancement of a-elements in elliptical galaxies 



7 




Fig. 3. The same of Fig. | but for the 6B-LD galaxy, the 
object with the lowest initial density and the longest du- 
ration of stellar activity. The total mass is Mt — 1O^M0, 
whereas the baryonic mass is Mt — lO^M©. The percent- 
age of Qf-enhanced stars is now 44% 



Fig. 5. Degree of enhancement in models with different 
shape and duration of the SFR. The dotted line is for 
B/A=0.1, the solid line for B/A=l, and the dashed line 
for B/A=10. See the text for more details 
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Fig. 4. Schematic SFRs of our simulations. See the text 
for more details. 



from which we may easily derive an analytical expres- 
sion for the SFR to be plugged into the chemical code. 
We have explored the following combinations of the pa- 
rameters B/A=0.1, 1.0, and 10 and ATsf=1, 2, 4, 6 
and 10 Gyr. The results are displayed in Fig. ||, which 
plots the present-day percentage of a-enhanced stars as a 
function of ATsf and the ratio B/A (the dotted line for 
B/A=0.1, the solid line for B/A=l, and the dashed line 
for B/A=10). It is soon evident that ATsf is the domi- 
nant parameter, whereas the ratio B/A plays a marginal 
role. With the prescription for Type la supernovae we 
have adopted, to fall below the border line of 50%, ATsf 
should be longer than about 6 Gyr and the stellar activ- 
ity should be uniformly diluted. This is indeed the case 
of 6B-LD model which finds confirmation by the present 
analysis. A counter-example is the case of the 7B model 
with Mt = 10^ Mq whose ATsf =^6.5 Gyr but the SFR is 
mostly confined between 2 to 4 Gyr with the peak at about 
3 Gyr. The reason for that is the initial density which is 
slightly higher than that of case 6B-LD (see Table |l| and 
[Chiosi fc Carraro (2001 ) for all details). The total du- 
ration of the SFR and not its shape drives the level of 
enhancement. 



build up in the galaxy. The galaxy under consideration has 
a total mass in stars of about 6.4 x 10^*^ M©. 

The following relationship holds between Mg and SFR 



B- 



Ms = 0.5 X A[l + y] ^TsF 

Pi. 



(4) 



5. Mean narrow band indices for galaxies 

In this section we evaluate the narrow band indices for 
our model galaxies in order to check how far they can go 
in reprod ucing important relat ionships such as the Mg2 
versus a ( Burstein et al. 1988| ) and the versus Mg2 
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Table 2. Real [Fe/H] for solar-scaled and a-enhanced According to Borges et al. (199E), the index Mg2 of a 
isochrones as a function of the initial chemical compo- 
sition [X,Y,Z]. 



"star" with enhanced content of Mg is given by the rela- 
tion 





r = 0.0 


r = 0.356 


Z Y X 


[Fe/H] 


[Fe/H] 


0.008 0.248 0.7440 
0.019 0.273 0.7080 
0.040 0.320 0.6400 
0.070 0.338 0.5430 


-0.3972 
0.0000 
0.3672 
0.6824 


-0.7529 
-0.3557 
0.0115 
0.3267 



lnMg2 = -9.037 5.795^^ 0.398 logg- 



T, 



off 



0.389 



"Fe' 




Fe' 


2 


"Mg" 




- 0.16 




-f 0.981 










Fe_ 



(9) 



and/or MgFe, see for instance Gonzalez (1993); [Trager 
et al. ( ^000a) Jb|). 

In presence of enhancement in a-elements one has to 



which holds for effective temperatures and gravities in the 
ranges 3800 < Teff < 6500 K and 0.7 < logg < 4.5. To 
make use of eq.(||), we need [Mg/Fe] as function of F. The 
correspondence between [Mg/Fe] and F can be obtained 
in the following way: from the pattern of a bund ances by 
Anders fc Grevesse (1989|) ; |G revesse (199l| ) and |Grevesse| 



and the iron content [Fc/H] by suitably dehnmg the en 
hancement factor F. For more details see Fantalo (1998) 
and [Tantalo fc Chiosi (2002[ ). 

Let us split the metallicity Z in the sum of two terms 



modify ] the relationship between the total metaliicity \ ^^^^^ (y^^^ ^^^^ approximated relation 



Z — ^ Xj + Xpe 



(5) 



where Icj are the abundances by mass of all heavy elements 
but Fe, and Xpo is the same for Fe. Recasting eq.(||) as 



Fo 



X: 
X^ 



Xf 



Xpc 



(6) 



and normalizing it to the solar values we get 



Z 



Xpc \ / X 



Xh / \Xpo/q vXh,o 



Xf 



Xfc 



(7) 



from which we finally obtain 



Fe 



Ip 



z 

z^ 



Ig 



X 

x;^ 



(8) 



which provides also the definition of F. Obviously for so- 
lar metallicity Z© and solar-scaled mixture (F = 0.0) 
eq.(||) yields [Fe/H]= 0.0. In general, when deahng with 
enhanced chemical mixtures, it suffices to determine F ac- 
cording to the degree of enhancement and simply re-scale 
the zero-point for the iron content of the Sun, i.e [Fe/H]=- 
F). Remarkably, enhancing a-elements simply reduces to 
decreasing [Fe/H]. 

Table |^ lists the initial chemical composition and the iron 
abundances corresponding to F=0 and F=0.356, the same 
as the stellar models and isochrones of Salasnich et al 
(2000|)]on which our indices are ultimately based. 

Once F is known for the stellar population under con- 
sideration, to proceed further we need a calibration ex- 
pressing the indices as a function o f it. To this purpose 
we make use of the calibration by Borges et al. (1995 ) 
for Mg2 in which the effect of enhancement is considered. 



0.8 



Mg 
Fe 



(10) 



According to which [Mg/Fe] =0.0 for F=0 (solar-scaled 
abundances), and [Mg/Fe]=0.444 for F=0.356. 



For all remaining indices we use the Worthey et al 
(1994D fitting functions in which the enhancement of 



a-elements is taken into account only via the re-scaled 
value of [Fe/H] of eq.(|). 

Looking at the response of indices to variations of the 
age, log(t) in Gyr, metaUicity log(Z/Z0), and degree of 
enhancement F of SSPs, we derive the following analytical 
relationships. 



li = (ai,iF + ao.i) log(t) 

+ (6i,,F + bo,,) log + co,,F + Si 



(11) 



where Ij for i=l to 4 stand for the Hp, Mg2, (Fe) and Mgb, 
res pectively. The coeffic ients are tabulated in Table ^ for 
the Borges et al. (1995 ) calibration. These analytical fits 
are accurate enough for the purposes of the present study. 

Given these premises, we derive the indices for our 
model galaxies as follows: 

(i) With the aid of the cumulative distributions of living 
stars per [0/Fe] bins (see the bottom panels of Figs. ||and 
^ and those for all remaining models not shown here for 
the sake of brevity) , we estimate the mean value of [O /Fe] 
weighed on the number of stars per bin and consider it as 
an estimate of the enhancement factor F. This is found to 
weakly correlate with the velocity dispersion of the galaxy 
F = 0.072 Zogcr-t- 0.134. 

(ii) From the history of star formation (top panels of the 
same Figures.) we estimate the mean age Tbs of the bulk 
stars. To do so one has to convert the rest-frame ages of 
the galaxy models into an absolute age by means of a cos- 
mological model of the Universe. To make things simple 
we adopt here the Freedman model with Ho=50 km s^^ 
Mpc~^, qo=0 and red-shift of galaxy formation Zfor = 5, 
which yield the maximum age for galaxies Tg=16.5 Gyr. 
The age Tbs is therefore Tbs = Tq — Tpsp, where Tpsp 
is the time at which the peak of star formation occurred. 
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bi,i 


bo,i 


CO,i 


3, 




-0.565 ± 0.134 


-1.356 ± 0.228 


-0.295 ± 0.140 


-0.673 ± 0.029 


0.205 ± 0.015 


2.89 


Mg2 


0.089 ± 0.034 


0.135 ± 0.056 


-0.072 ± 0.027 


0.165 ± 0.030 


-0.093 ± 0.010 


0.10 


(Fe> 


0.488 ± 0.132 


0.845 ± 0.580 


0.357 ± 0.113 


1.725 ± 0.185 


0.898 ± 0.026 


2.18 


Mgb 


1.103 ± 0.234 


1.676 ± 0.675 


0.115 ± 0.298 


2.011 ± 0.359 


0.806 ± 0.066 


2.09 



{a = 0,lB6 b=-0.ia5) 

= 0.092 b==0.027; 
(a = 0.062 b = 0.004) 

3x<Z> (a-0,092 b-0-066) 

Composite (a = 149 b = -0.069 



Data 
(Z) 




A A - ■' _A 



Fig. 6. The Mg2 versus velocity dispersion as- The small 



open circles are the observational data from Trager et al 
(2000a[^) and the solid line is their least squares fit. The 
open and filled triangles are the theoretical models of type 
B listed in Table ^ for the maximum and mean metallic- 
ity, respectively. Their least square fits are shown by the 
dashed (ZMax) and dotted ((Z)) Hues. The three empty 
stars show the same, however limited to the three most 
massive galaxies, but for a metallicity arbitrarily increased 
to 3 X (Z). The thick dashed line shows the composite case 
in which the metallicity of the three most massive galax- 
ies is 3 X (Z) and that of the less massive ones is (Z). The 
filled squares (ZMax) and boxed-filled triangles (langleZ) 
show the case of models A. 



These quantities are listed in columns (5) and (6) of 
Table g 

(iii) Finally, we make use of relations (^l]) to derive the 
indices H^, Mg2, (Fe), and Mgb for both the maximum 
and mean metallicity (age Tbs and F are the same). The 
results are listed in Table ^. 

The Mg2 vs a relation. This is shown in Fig. H an d 



compared with the observational data by Trager et al 
(20004|) for a sample of local EGs. The least squares 
fit of the observational data yields the relation Mg2 — 
0.186 logo- — 0.125 (where a is in kms~^). The theoretical 
results are plotted separately for each of the metallici- 



ties in use: the filled triangles are for (Z) (a more realis- 
tic approximation of reality) and the open triangles are 
for Ziviax- The least squares fit of the data yields Mg2 = 
0.063 logo- -HO.OOl for (Z) and Mg2 = 0.094 logcr-H0.063 for 
ZMax- The slope is about half of what observed. Where is 
the source of disagreement? There are several candidates: 
(i) the velocity dispersion, (ii) the calibration in use, (iii) 
the set of parameters Tbs, Z and F assigned to the model 
galaxies. Limiting the discussion to the range of massive 
models for which data exist, we note: 

(i) The velocity dispersion is not likely because it should 
be decreased by a factor of two and in any case the relation 
would simply horizontally shift without changing its slope. 

(ii) The calibration is more uncertain, even though the 
one in use here allowing for the maximum effect. In this 
context one should, however, keep in mind that, contrary 
to what one may guess, increasing F (higher enhancement) 
does not imply a stronger Mg2 because of the counter- 
reaction by [Fe/H]. 

(iii) Assuming that our velocity dispersions and calibra- 
tion are correct, inspection of the weight of Tbs, T, and 
Z on building up the final result reveals that the metal- 
licity dominates. The three empty stars show the effect of 
arbitrarily increasing (Z) by a factor of 3, limited to the 
three most massive galaxies. This implies that the Mg2 vs 
a relation reflects the mass-metallicity relation, rather a 
sequence in the enhancement level. 

Are the mean (and maximum) metallicity reasonable? 
Is there any plausible effect making a massive galaxy more 
efficient in building metals? With the IMF adopted here 
and in the dynamical models, the production of metals 
is scarcely efficient because of the low Increasing this 
for all galaxies would simply shift vertically the Mg2 vs 
a relation. However, if C increases with the galaxy mass 
(top-heavy IMF), theory and data could be reconciled. 
Although this possibility has often been advanced by sev- 
eral authors in different astrophysical contexts - see the 
arguments given by Chiosi et al. (199^ ) and Larson (1998[ ) 
- further exploring the subject is beyond the scope of this 
paper. Suffice to conclude that dynamical models almost 
match the observational Mg2 vs a relation provided some 
fine tuning of the chemical yields per stellar generation 
is applied. Just for the sake of curiosity, we derive the 
relationship one would obtain by letting the mean metal- 
licity to increase from the present values to 3 x (Z) as the 
galaxy mass mass increases from 10^ to 5 x 10^^ Mq. We 
get Mg2 = 0.149 logcr — 0.070, which is about right. 

The Hfj versus Mg2 and MgFe planes. The index-index 
planes are customarily used to infer age and metallicity of 
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Table 4. Theoretical indices for all the models under considerations. The calibration in use is by Borges et al. (1995). 
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lA 


1 X 10" 


0.020 


0.005 
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16.3 
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0.26 
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3.09 


1.79 


0.14 
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0.002 
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1.33 
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2.33 


2.99 


2B 
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0.15 
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1 X 10^ 
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0.01 
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0.04 
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1.39 
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1 X lO'' 
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0.12 


10.0 


2.01 


0.12 
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2.14 


0.08 
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1.78 


6B-HD 


1 X 10^ 


0.006 


0.002 


0.19 


12.0 


1.81 


0.15 


1.99 
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Fig. 7. The Hf3 versus Mg2 and MgFe planes. The mean- 
ing of the symbols is the same as in Fig. 0. We have also 



added the small sample of EGs in Fornax by Kuntschner 
(2000|)](big filled circles). The line is the linear fit of all 
observational data. 



early-type galaxies ( 


Gonzalez 1993; 


Bressan ct al. 1996; 


Trager et al. 2000a 


b|; Kuntschner 2000|, 2001). In Fig. |7| 



we present the vs Mg2 and vs MgFe planes left 
and right panels respectively, for our models, single stel- 
lar populations of similar metallicity, and observational 
data for a small s ample of EGs in the Fornax cluster by 
Kuntschner (2000 ) and the wider sample of nearby galax- 
ies by Trager et al. (2000a ) . The following remarks can be 
made: 

(i) Our models lie on a strip which roughly correspond to 
that covered by SSP of different metallicity and age in the 
range say 15 to 8-10 Gyr. This expected because all models 
have terminated their star formation far in the past. Even 



the most extreme case has ceased star formation 6 Gyr 
ago. 

(ii) Also in these planes we notice the same kind of dif- 
ficulty encountered with the Mg2 vs a relation, i.e. our 
model galaxies do not make enough metals so that even 
the most massive ones only marginally reach the region 
crowed by the data. The same kind of arguments brought 
about to cure the Mg2 vs a relation apply here. Increasing 
the metal content of the model galaxies is less of a prob- 
lem. 



6. Concluding remarks 

In this paper we have derived the present-day level of 
enhancement in a-elements in the stellar content of the 



N-Bod y Tree-SPH models of EGs by [Chiosi fc Carraro 
|(200lD , cast light on the main cause of it, derived the mean 
spectral indices of the model galaxies, and compared them 
to the observations with particular attention to the Mg2-cr 
relation and the - Mg2-[MgFe] planes. 

Combining the results of Chiosi fc Carraro (2001| ) sum- 
marized in Section 3 and those obtained here the following 
remarks can be made: 

(i) The star formation history of EGs depends on their 
total mass and initial density, i.e. the physical conditions 
of the environment in which they have formed. High- 
density systems, owing to their very early and short star 
forming activity would now appear as old galaxies, whose 
stellar content is enhanced in a-elements. Low-density sys- 
tems of high mass would resemble the previous ones, their 
stellar content should, however, span a broader range of 
ages. Most of the component stars should be enhanced in 
a-elements but with a broader distribution of abundance 
ratios [a/Fe]. At decreasing mass the star formation activ- 
ity should become more diluted in time, the stellar popula- 
tions should span wider and wider age ranges, and the level 
of enhancement should progressively tend to solar. The 
extreme case is reached with the low-mass, low-density 
systems which could undergo many episodes of star for- 
mation up to very recent times, thus resembling young 
objects with solar or sub-solar partitions of elements. It is 
worth recalling here that F or equivalently [O /Fe] is found 
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to weakly increase with the velocity dispersion (mass) of 
the galaxy. 

(ii) The scenario sketched above is perhaps confirmed 
by the many observational studies based on different kinds 
of diagnostic. For instance Kmitschner fc Davies (199^ ); 
Kuntschner (200C, 2001) studying early-type galaxies in 
the Fornax cluster find the EGs to be coeval and to have 
metallicities varying from solar to three times solar, and 
the ratio [Mg/Fe] -inferred from the indices Mg2 and Fe3- 
to weakly correlate with the velocity dispersion. In con- 
trast the data for nearby fi eld galaxies by [ Gonzalez (1993| ) ; 
Bressan et al. (199^ ), and Trager et al. (2000a ,|h|), seem to 
indicate a much broad range of ages as perhaps compati- 
ble with our model galaxies of low initial density. Finally 



we call attention on the recent result by Poggianti et al 
(2 00l|)]for dwarf galaxie s con firming the suggestio n made 
by [Bressan et al. (1996| ) and [Trager et al. (20004 H) of a 
recent star formation in these systems. 

(iii) The Mg2 vs u relation mainly reflects the 
metallicity-cr-mass sequence rather than enhancement in 
a-elements and/or age. 

(iv) Although the dynamical models and companion 
chemical models we have discussed do not fully match the 
high metallicities inferred from other diagnostics and some 
improvement is required, they however seem to offer a 
coherent scenario accounting for many observed properties 
of EGs, thus perhaps indicating that we are on the right 
track. 
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